High temperature expansion applied to fermions near Feshbach resonance 
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We show that, apart from a difference in scale, all of the surprising recently observed properties 
of a degenerate Fermi gas near a Feshbach resonance persist in the high temperature Boltzmann 
regime. In this regime, the Feshbach resonance is unshifted. By sweeping across the resonance, 
a thermal distribution of bound states (molecules) can be reversibly generated. Throughout this 
process, the interaction energy is negative and continuous. We also show that this behavior must 
persist at lower temperatures unless there is a phase transition as the temperature is lowered. We 
rigorously demonstrate universal behavior near the resonance. 

PACS numbers: 03.75.-b,03.75.Ss,34.50.-s 



At present, much experimental activity is concentrated 
on degenerate quantum gases near a Feshbach resonance 
[H [3- C3- G3- 01 where the nominally weak effective inter- 
actions are strongly enhanced. By applying a magnetic 
field, which moves the energy of a bound molecular state 
relative to the scattering continuum, the interactions can 
be tuned by by many orders of magnitude. The scatter- 
ing length diverges to negative/positive infinity when the 
molecular energy is infinitesmally above/below threshold. 
At resonance, the scattering cross-section is limited only 
by unitarity (er = Air/k 2 , where k is the relative momen- 
tum of the two atoms) and is universal (independent of 
the microscopic physics). The universal nature of the uni- 
tarity limit leads to a challenging many-body problem, 
since there are no small parameters readily identifiable 
for the application of perturbation theory. 

Recently, a sequence of beautiful experiments on Fermi 
gases have found unexpected and dramatic behavior near 
unitarity; including a universal interaction energy 0, Q 
and the reversible interconversion of atoms and molecules 
Q . The former demonstrates the effects of unitary scat- 
tering on bulk properties, and the latter suggests the pos- 
sibility of an equilibrium phase with atoms and molecules 
in chemical equilibrium. Several of the experimental re- 
sults are not well understood, and some even appear to 
be mutually contradictory Q- In the following, we shall 
summarize these experiments and the fundamental ques- 
tions which they raise. We refer to the range of magnetic 
fields on cither side of the resonance where the scatter- 
ing length a sc is positive or negative as the "+'ve" and 
"— 've" side of the resonance. 

(I) Universal interaction energy: Thomas's group at 
Duke 1] has studied a gas containing two spin states 
of fermionic 6 Li (which has a Feshbach resonance around 
855 Gauss). The experiment is performed on the nega- 
tive side at 910 Gauss. Upon release from the trap, the 
Fermi gas undergoes anisotropic expansion for tempera- 
tures between 0.1 and 3.5 Tp. This anisotropy can be 
explained by collisional hydrodynamics with a universal 
interaction energy proportional to the Fermi energy £p. 



(II) Properties near resonance: More recently, Salomon's 
group directly measured the interaction energy of the 
same 6 Li system, at temperatures between 0.5 and 1 Tp. 
(a) Crossing the resonance from the negative side, they 
find that the interaction energy e- ln t remains negative and 
continuous across the resonance, despite the expected in- 
finite jump of the scattering length. They find an inter- 
action energy similar to that of the Duke group, convinc- 
ingly demonstrating that ej n t remains roughly constant 
over the temperature range O.ITf to Tp. (b) Approach- 
ing the resonance from the opposite side, ei n t is positive 
but drops to a negative value at a field of 700G; roughly 
150G before the resonance is reached, (c) The three- 
body recombination rate is maximal at 700 Gauss, rather 
than at 855 G; a result consistent with the work of Ket- 
terle's group at MIT |3j. (d) The anisotropy of the ex- 
pansion (which is a measure of the interaction strength), 
is maximal at 855 G. While (a) through (c) seem to indi- 
cate that the resonance is shifted from 855 G to 700 G by 
many-body effects; (d) is consistent with an un-shifted 
resonance. 

(III) Conversion between atoms and molecules: Very re- 
cently, Jin's group at JILA has studied a gas consisting 
of two spin states of fermionic 40 K £| . They show that 
molecules are produced as one crosses the resonance from 
the — 've to the -f-'ve side. Moreover, this process is re- 
versible. This experiment, performed at T = O.ITf finds 
no significant shift of the original resonance. 

Here, we focus on the Fermi gas near a Feshbach res- 
onance. Far from resonance, the interaction energy is 
£int = 5^t n l' where g = Ann 2 a sc / 'M , nj and are the 
number densities of the two spin components. (We shall 
consider the case nj = = n/2.) A key question is how 
this non-universal form of e- m t turns into a universal func- 
tion near resonance, (see (I) and (II))- Moreover, how 
does this function reflect the existence of molecules and 
what signatures do these molecules have if they exist? 
Equally important is whether the original resonance is 
shifted in a many-body medium, (see (II))- To use meth- 
ods where errors can be estimated precisely, we consider 
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the Boltzmann regime. In this case, physical quantities 
can be calculated systematically through a high tempera- 
ture series expansion; and yet the issues of the emergence 
of universal behaviors at unitarity still remain. 

One might wonder whether the phenomena in the 
Boltzmann regime have any relevance to current ultralow 
temperature experiments. The connection is simply that 
in the absence of any phase transition in the normal state, 
the thermodynamic functions in the Boltzmann regime 
are analytically continuations (in T and n) of those in the 
degenerate regime. This analytic continuation imposes 
strong constraints on the phase diagram, allowing one to 
qualitatively understand the physics at all temperatures. 
In fact, apart from a difference in scale, the exact re- 
sults in Boltzmann regime show all of the features discov- 
ered in experiments (I) , (II) , and (III) . Surprisingly, the 
high temperature results, which should work well when 
n T A 3 = n x A 3 = n\ 3 /2 < 1 (where A = h/y/2-KMk B T 
is the thermal wavelength and ks is Boltzmann's con- 
stant,) agrees reasonably well with the ENS experiments 
@ performed at nA 3 /2 = (4/30F) (7>/T) 3 / 2 - 1.6 (see 
point (E), later). In the following, we derive the interac- 
tion energy density in the Boltzmann regime, then draw 
a series of conclusions labeled below as (A) to (D). 

Energy density in the Boltzmann regime: At 
high temperatures, or low densities, the grand parti- 
tion function Z = r £re~( H ~ ,J ' N ^ kBT can be expanded 
in the fugacity z = e^/ kBT 0, 0; and to the sec- 
ond order in z, where interaction effects first appear, 
is Z = + 2\f2 (Vz 2 /\ 3 ) b 2 , where the superscript 
"o" denotes quantities for non-interacting systems, and 



cient, 



e 1 ) is the second virial coeffi- 



b 2 = Y / e lEbl/kBT 



dk d8t(k) 
7r dk 



(1) 



where *yt = 21 + 1, the sum is over all integers £ = 0, 1, .., 
Eb is the energy of the two body bound state, and 8g(k) 
is the phase shift of the £-th partial wave|g. 

A standard thermodynamic calculation [9| shows that, 
to lowest nontrivial order in nA 3 , the energy density is 
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where e^in and e; n t are respectively the kinetic and inter- 
action energy densities; 
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Since the contributions of the partial waves with £ > 1 in 
eq. are smaller than those of the s-wave by a factor 
nA 3 , we shall ignore them in our subsequent discussions. 

Far from resonance the phase shift is S(k) = —a sc k for 
small k. Due to the Gaussian cutoff in eq. QJ, no errors 



are made by using this expression for all fc. In the absence 
of bound states, we then have b 2 — —a sc /(y2X), and 
hence Tdb 2 /dT = b 2 /2. Equation then reduces to 
the usual expression eint = gn^n^, with n<\ = n± = n/2. 

Energy levels and interaction energy : To eval- 
uate eint near Feshbach resonance we use the expression 
for the phase shift valid for fc << where b is the 
range of the potential |Toj 



fccot <5(fc) 



r a k 2 



(4) 



where r (~ b) is the effective range and a sc is the s-wave 
scattering length. As a function of magnetic field, a sc 

varies as a sc = at, g (l — B A B B J , where B and AB are 

the location and the width of the resonance, and a^g is 
the background scattering length away from resonance. 
The region B > Bo and B < Bq where a sc < and 
> are referred to as the "— 've" and "+'ve" side of 
the resonance. For temperature range such that thermal 
wavelength is larger than the range of the potential, A > 
b, we can use eq.Q in eq.Q) where the integral is cutoff 
at A -1 . Near resonance, X/a SCl r /a sc —>■ 0, and eqs.Q 
and 10} imply that up to terms of order r Q / a sc , 



b 2 = Y, e |£ " l/fcBT - (1 - erf (x)) e* 



(5) 



where x = A/(v27ra sc ), and erf(x) is the error function. 
We can then evaluate ej„t using eq.® and J^J. At the 
same time, it is useful to look at the effect of phase shifts 
on the energy levels. In a box of size R, the wavevector 
fc is changed from its non- interacting value fc = £tt/R 
(£ = 1, 2, 3..) to fc ~ fc — 8{k )/R through the boundary 
condition sin(fci? + 5(k)) = at a large distance R. Scat- 
tering states and bound states correspond to real and 
imaginary fc solution of this boundary condition. The 
energies of lowest few states in the center of mass frame 
(E = h 2 k 2 /M) are shown in fig. 1. As one passes through 
the resonance from the — 've side, the lowest state in the 
continuum turns into a bound state, causing S(k = 0) to 
change abruptly from to 7r. The corresponding inter- 
action energy e m t calculated from eq.JSJ) and eq.® are 
shown in fig. 2. Despite the simplicity of the calculation, 
considerable information can be deduced: 

(A) Approaching the resonance from the — 've side; 
follows a "negative branch" e[ nt ^ which is negative and 
decreases monotonically. (See fig. 2). e[ nt ' ) evolves from 
the (temperature independent) non-universal form gn\n^ 
far from resonance to a (temperature dependent) univer- 
sal form — eo = — (3rtfcsT/2)(nA 3 /2 3 / 2 ) at resonance, and 
continues on to the +'ve side. This universal form, which 
follows from the fact that b 2 = 1/2 and db 2 /dT = at 
resonance (x = 0)0, is the high temperature analog of 
universal interaction found in (I) and (II). 

(B) Despite the change in sign of the scattering length, 
the interaction energy e int remains negative across the 
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FIG. 1: (Color Online): Energy levels of a two body system in 
the center of mass frame calculated from the phase shift and 
the boundary condition in the text, with R — 14/im, Bo = 
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FIG. 2: (Color Online): Interaction energy. The ener- 
855G, AB = 325G, and a bg = -120nm. The vertical scale is S ies 4tt increase monotonically as the field is increased, 
expanded to isolate individual energy levels. The dotted lines reaching the universal value ±e at resonance, where e = 
are energy levels of a non-interacting system, E ( ° ] = E n 2 , (3nk B T/2)(n\ 3 /2 3/2 ). The negative branch, e[J, is continu- 



E = h 2 (n/R) 2 /M, n = 1,2, 3. Passing from B > B to B < 
Bo, the lowest state in the continuum becomes bound. Light 
(yellow) upright and dark (red) inverted triangles indicate 
states at the same field with and without the bound state 
occupied. The interaction energies due to thermal occupation 
of these states, and those marked by circles and squares, are 
shown in fig.2 by the same symbols. The inset shows the 
behavior of the scattering length. 



resonance (shown in fig. 2) because of the thermal popu- 
lation (rib) of bound states which exists when a sc is pos- 
itive, 



ous across the resonance, while ett^ will jump to ek* if the 
bound state is occupied (say, at the field labeled by the trian- 
gle). Inset shows the temperature dependence: solid, dashed 
and dotted lines correspond to T = 1, 10, 100/iK. 
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Although proportional to the small factor n\ 3 , (rib) is 
macroscopic in the thermodynamic limit. Equation © 
in turn implies that in a bulk system, the original reso- 
nance (at field Bo) can not be shifted to the positive side 
(to B\ < Bo) at a lower temperature unless there is a 
phase transition for all magnetic fields between B\ and 
Bq where (rib) disappears as temperature is lowered from 
the Boltzman regime. (See fig. 3). So far, such a phase 
transition has not been observed. Should future exper- 
iments rule out such a transition in 6 Li, one must then 
conclude that the resonance in the ENS experiment [2j is 
not shifted. The absence of a shift of the resonance is 
also consistent with the findings in ref. T] for 40 K. 
(C) Approaching the resonance from the +'ve side; if the 
bound states are not occupied, £i n t will follow a "positive 
branch" ej^ > which increases monotonically, evolving 
from gn^ni to e = (3nfcsT/2)(nA 3 /2 3 / 2 ) at resonance 
(see fig. 2). However, a gas containing only scattering 
states is not in true equilibrium, as the latter requires 
that all states, including bound states, are thermally pop- 
ulated. As pointed out by Petrov the three-body col- 
lisions which convert scattering states into bound states 



FIG. 3: Schematic phase diagrams. Because the original 
resonance (-Bo) is un-shifted at high temperature, a shift in 
resonance at low temperature (-Bi < Bq) will imply a phase 
boundary as shown in (a), which means that (rib) will disap- 
pear (or appear) as temperature is lowered (or raised) along 
the path indicated by the double-headed arrow. The phase 
diagram for an unshifted resonance is shown in (b). 



can only lead to chemical equilibrium if the released en- 
ergy is insufficient to eject the particles from the trap. 
If one is sufficiently far from the resonance, then the 
molecular binding energy is larger than the trap depth, 
and the three-body collisions instead lead to loss. The 
sudden drop in the interaction energy near 700G in the 
ENS experiment 2j is consistent with the production of 
trapped molecules. Moreover, Jochim et al. have di- 
rectly observed the equilibration of atom/molecule pop- 
ulations near this field. This scenario also predicts that 
the development of an equilibrium molecular population 
will coincide with a peak in the three body loss rate. This 
peak is observed at ENS @, MIT 0, and Innsbruck 0. 
(D) The extension of e[~^ to the positive side of the 
resonance is due to the population of the bound state 
which becomes available there. If the system is brought 
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across the resonance adiabatically, quasi equilibrium is 
maintained and the process is reversible. Consequently, 
molecules generated on the positive side will turn back to 
atoms when the system is brought back to the negative 
side. This is consistent with the experiments on 40 K 
(E) The experiment at ENS was performed at tempera- 
ture T — 3.5/iK and degeneracy factor T/Tp — 0.6, cor- 
responding to nA 3 /2 = 1.6. Using the same temperature, 
we have plotted the ratio eint/ekin in fig. 4 for T/Tp — 
1.2, 0.6, and 0.4, corresponding to n\ 3 /2 = 0.6, 1.6,3.0. 
Although our result on interaction energy should only be 
accurate when nA 3 /2 < 1, we also plot it at higher phase 
space densities (i.e. extending it to regions nA 3 /2 > 1) 
to indicate its temperature dependence. Moreover, these 
extensions are the leading term of the interaction energy 
in powers of nA 3 /2. On the positive side all three of our 
curves are consistent with the experimental data, while 
on the negative side the T/Tp — 1.2 curve is closest. The 
fact that the B > Bq data matches this higher tempera- 
ture curve may be due to higher order terms in the high 
temperature expansion or systematic differences in the 
density/temperature of the sample on the two sides of 
the resonance. In any case, it is clear is that a Fermi 
gas in the Boltzmann regime exhibits all the phenomena 
seen in current experiments. Equally important is the 
fact that the exact high temperature results near unitar- 
ity forces one to conclude the existence of a phase where 
atoms and molecules are in chemical equilibrium. 

Having shown that the behavior of the quantum gas 
in Boltzmann regime has all the characteristics seen in 
experiments at lower temperatures, we reiterate that the 
two regimes are connected by the analyticity of the ther- 
modynamic functions. Despite the divergent scattering 
length, at small fugacity z (or equivalently small den- 
sity nA 3 ) the probability of particle collisions is small, 
leading to our systematic perturbative scheme. 
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FIG. 4: (Color Online): The ratio eintAw for T = 3.5^K. 
Squares (circles) show data from ref. [3 which agrees with 
e + , (e - ). Triangles do not fall on either curve and probably 
reflect a nonequilibrium situation. The dashed, solid, and 
dotted lines are for T/Tf = 1.2, 0.6, 0.4 respectively. 



